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Abstract 

We study the deterministic dynamics of a periodically driven particle in the under- 
damped case in a spatially symmetric periodic potential. The system is subjected to 
a space-dependent friction coefficient, which is similarly periodic as the potential but 
with a phase difference. We observe that frictional inhomogeneity in a symmetric 
periodic potential mimics most of the qualitative features of deterministic dynamics 
in a homogeneous system with an asymmetric periodic potential. We point out the 
need of averaging over the initial phase of the external drive at small frictional in- 
homogeneity parameter values or analogously low potential asymmetry regimes in 
obtaining ratchet current. We also show that at low amplitudes of the drive, where 
ratchet current is not possible in the deterministic case, noise plays a significant 
role in realizing ratchet current. 
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1 Introduction 

Directed particle transport in asymmetric periodic (ratchetlike) potentials 
without the application of any obvious external bias but aided by thermal 
fluctuations has been widely studied both in overdamped and underdamped 
conditions [1 ,2,3,4.5] . This effect, called the ratchet effect, has also been ob- 
tained in underdamped symmetric periodic (non-ratchetlike) potential sys- 
tems, by considering either the temperature [6II71I8I9] or the friction coefficient [ToyTT] 
to be space dependent. Ratchet current can also result in a symmetric peri- 
odic potential when the system is driven by an externally applied zero-mean 

* Corresponding Author 

Email address: mangal Onehu . ac . in (Mangal C. Mahato 2 ). 



Preprint submitted to Elsevier 



6 April 2010 



temporally asymmetric forcing aided by Gaussian white noise [12] or by the 
explicit application of nonthermal noise |13(Ti] . The ratchet effect is being 
studied intensively for about two decades and its applications envisaged in 
periodic systems that are common in many areas of natural sciences [To]. 

All the above ratchet models use the presence of (Gaussian or colored) noise 
(fluctuations) to obtain directed transport. In these models noise is essential. 
However, the study of deterministic ratchets, pioneered by Ref.[16], has also 
contributed significantly to the understanding of the subject. These determin- 
istic ratchets, unaided by noise, are shown to yield current in overdamped [17)115] . 
underdamped [TlfTo1fT9T2(^ as well as in Hamiltonian[29|[30||3T] 

periodic potential systems, and also in overdamped quenched disordered [32|M] 
systems. In these systems net current results, without the presence of applied 
nonzero average forcing or asymmetric fluctuations, due to the presence of 
various regular transporting or chaotic attractors depending on the initial con- 
ditions for given system parameter values. In Ref. [H] the symmetry criteria 
for the realization of ratchet current have been discussed in detail. 

Owing to the complexity of the dynamics, for example the simultaneous pres- 
ence of periodic and chaotic attractors, there is a need for proper choice 
of initial conditions and ensemble averaging over them to obtain realistic 
averages [TUl 1 9 II 2 9] . This fact has been exemplified in an (zero mean) ac mod- 
ulated periodic Hamiltonian system [31] . It has been shown there that slower 
modulations of the symmetric periodic potential lead to asymmetric access 
(in momentum) of chaotic regimes in the phase space resulting in a giant net 
particle transport in one direction. The direction and magnitude of the current 
depends on the initial phase of the modulating ac drive. Of course, the net 
current must disappear if averaged over the initial phase. 

The driven deterministic systems show frequent ensemble- averaged current 
reversals as the drive parameters are changed. Mateos identified bifurcation 
from a chaotic to a periodic regime as the mechanism for the average- current 
reversals |20j in these systems. However, the change in direction of individual 
single particle trajectories could be related to phase locking phenomena [21] 
due to the presence of various velocity attractors [21]. The average current di- 
rections are thus sensitive to initial conditions [38] which allow to explore these 
attractors selectively in the phase space. In a recent further investigation [19] 
it has been conjectured that, in general, bifurcations from chaos to periodic 
behaviour could be associated with abrupt changes (and not necessarily re- 
versals) in the ensemble averaged currents. 

In Ref . [31] , the effect of different control parameters were analysed in the iner- 
tial limit. It was found that the control mechanisms were associated with the 
fractal nature of the basins of attraction of the mean velocity attractors and 
that small perturbations of the control variable could produce drift reversals. 
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However, the effect of the presence of a weak subharmonic component in the 
ac drive field on the phase locked dynamics of a ratchet is to suppress chaos 
and stabilize regular orbits over large range of driving amplitudes |3 7]. 

The transport properties of inertial systems depend sensitively on the (con- 
stant) friction parameter [26. 27. 28J. The addition of noise to the system in 
the low damping regime makes the system dynamics robust against initial 
conditions [27]. Moreover, in a well defined damping window the efficiency of 
ratchet current generation becomes appreciable even for small applied peri- 
odic field amplitudes. In Ref. [2B] the effect of noise on individual transporting 
deterministic trajectories is analysed. A small noise seem to have little effect 
on some trajectories in one direction whereas it destroys trajectories in the 
opposite direction resulting in a robust finite ratchet current. 

In most of the above works (barring a few, for instance[14,29,31j), particle 
motion in periodic but asymmetric potentials is considered; for example a 
system with potential V(x) = Vo(cosa; + 6cos2x) driven by a periodic field, 
F(t) = acos(ut + 4>o) in uniform friction media. In these studies the po- 
tential asymmetry (b ^ 0) is the primary cause of net particle drift. In 
the present work, we consider a symmetric periodic potential, for example 
V(x) = — Vo sinx, but in the presence of a space dependent friction coefficient 
j(x) = 7o(l — A sin(x + 0)) (as in [10]). A massive charged particle moving in 
a periodic ionic lattice and in a medium with material density (friction) pro- 
file created by a stationary pressure wave can be thought of as a mechanical 
illustration of such systems [T0~|ITTj . However, such systems have been studied 
earlier in a resistively and capacitatively shunted junction circuit model of 
small Josephson junctions (3900]. The frictional inhomogeneity term emulat- 
ing there the " cos 9" term representing the coupling between the quasiparticle 
tunneling and Cooper pair tunneling across the Josephson junction. In Ref. 
[59] the phase space trajectories of this system are presented in detail. 

In the present work we show that all the features of deterministic inertia 
ratchets can be observed even in the symmetric periodic potential system in 
an inhomogeneous medium. The frictional nonuniformity of the medium with 
phase shift <f>(^ 0, n) essentially emulates the effect of the asymmetry of the 
potential considered in earlier cases. Moreover, we show in our present inhomo- 
geneous dissipative system that, for small inhomogeneity (roughly A < 0.6), in 
order to get sensible ratchet current one also needs to average over the initial 
phase 0o of the drive as in the Hamiltonian system [291I3T] . Though ratchet 
current is a steady state (t — > oo) phenomenon one is likely to obtain a finite 
ratchet current, an obvious erroneous result, even at zero potential asymme- 
try in a uniform medium without this averaging. Interestingly, the bifurcation 
diagrams with and without the frictional inhomogeneity are qualitatively in- 
distinguishable. However, the stroboscopic Poincare plots differ qualitatively 
in the driven system as the inhomogeneity parameter A is varied. Links Vi 
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The presence of noise (which is ubiquitous in natural environments) smooths 
out the current fluctuations so numerously present in the deterministic case 
[2"5f4"Tj . In Ref. [H] a matrix continued fraction method is presented to ob- 
tain ratchet current in noisy periodically driven inertial systems. We further 
show in our inhomogeneous system that at low drive amplitudes noise plays 
a crucial role in getting ratchet (mean) current, that is conspicuously absent 
in the deterministic drive case, though at larger drive amplitudes one obtains 
ratchet current without the presence of noise. 



2 The model 



As mentioned earlier, in this work, we consider the motion of a particle in 
a periodic potential V(x) = — Vosin(fcr) which is symmetric in space (about 
kx = (2n + l)7r/2, n = 0, ±1, ±2, ...). The friction coefficient 7(2;) = 7o(l — 
\sm(kx + 0)) is periodic with the same periodicity as the potential but has 
a phase difference. The system is driven periodically by an external periodic 
forcing F(t)=a cos(ut + <f> ). We study the system in the deterministic regime 
i.e. in the absence of noise and also at the end in the presence of noise. 

The one dimensional equation of motion of a particle of mass m is given by 
the Langevin equation, 

The corresponding equation of motion in the presence of noise is given by 
d 2 x .dx dV(x) 



m— - = —'yyx 



dt 2 dt dx 



+ F(t) + yfi(dft(t). (2) 



In Eq.2, T is the temperature in units of the Boltzmann constant ks- The 
Gaussian distributed fluctuating forces £(t) satisfy the statistics: < £(£) >= 0, 
and < ) >= 26{t - 1 ). For convenience, we write down Eq.l and Eq.2 

in dimensionless units by setting m — 1, Vo = 1, k — 1. In Eq.2, T = 2 
then corresponds to an energy equivalent equal to the potential barrier height 
at a = 0. In terms of the reduced variables denoted again now by the same 
symbols, the two equations can be written as 

d 2 x dx 

= -7(^ + cosa; + F(t), (3) 



di 'jo , , doc 



and 

dx 1 

d1> " 7(X) dt +cosx + F ^ + V7C*0 T £(*)> ( 4 ) 
where 7(2;) = 7o(l — A sin(x + 0)). Thus the periodicity of the potential V(x) 
and also the friction coefficient 7 is 2n . The potential barrier between any 
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two consecutive wells of V(x) persists for all field amplitudes a, (0 < a < 1) 
and it just disappears at the critical field value a = a c = 1. Notice that for 
= rnr, n = 0, ±1, ±2, .... the deterministic Eq.3 should yield exactly the same 
trajectory but in the opposite direction on letting x — > x — n/2 and t —¥ —t 
and in the presence of a thermal noise Eq.4 should give the same ensemble 
averaged current but again in the opposite direction [14]. Thus in order that 
a ratchet current appears the symmetry of the system must be broken using 
cf) 7^ mi. 

The noise variable, £ in Eq.4, satisfies exactly similar statistics as earlier. 
Eq.3 and Eq.4 are solved numerically using the Heun's method for solving 
differential equations. The particle trajectories are obtained for given initial 
conditions for fixed values of a and 70. Also, the steady state mean velocity of 
the particles are calculated. 

For the deterministic case the average velocity for a single particle (with a 
given initial condition) is obtained by the relation [24] 



where t max is the maximum time for which the particle is allowed to evolve 
and tt r an is the duration of initial transients that we remove. 

Also the mean velocity averaged over an ensemble of N particles is defined as 



For the particle trajectory in the presence of noise, the steady state mean 
velocity v of the particle is obtained as 



where the average (...) is evaluated over a large number of trajectories. 



3 Numerical Results 

The phase space trajectories of the deterministic Eq.3 was investigated earlier 
in Ref. [39]. Also, for symmetric systems (A = 0, in Eq.3) the regimes of 
chaotic (nontransporting) and periodic motion were given in detail in Ref. 
[42] . In the present work the particle is allowed to evolve for a long time 
so that the system reaches a steady state condition to investigate ratchet 
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Fig. 1. Plot of particle trajectory for a = 0.5 and r = 10 (main figure) and the 
corresponding phase space trajectory; x(0) = ir/2, v(0) = 0, 70 = 0.12, (ft = 0.35, 
A = 0.9. 
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Fig. 2. Variation of mean velocity (v) with amplitude a for a single particle; v(0) = 0, 
x(0) = tt/2, 70 = 0.12, <ft = 0.35, A = 0.9 . 

effect. We solve Eq.3 (for the deterministic case) and Eq.4 (when noise is 
included) using the second order Heun's method with time step size At = 
0.001. For the calculation of deterministic average velocity (Eq.5) we take 
tmax = 10 6 and ttran ~ 10 4 . Throughout our work, unless otherwise stated, we 
keep A = 0.9, 70 = 0.12 and (ft = 0.35. We partly repeat the investigations done 
earlier [20|21f24] replacing the asymmetric periodic potentials by a symmetric 
periodic potential but instead of a uniform-friction-medium we consider here 
a periodically varying friction coefficient and present our results below. 

3.1 Phase locked dynamics and sensitivity to initial conditions 

The dynamics of the particle is critically dependent on the amplitude a and 
frequency u of the zero-mean external periodic drive [31]. The dissipative sys- 
tem goes to the only fixed point attractor at v = 0, x = (2n — l)ir/2 when the 
system is not driven (a = 0). For a given (non-zero) amplitude and frequency 
of the drive the individual particle trajectories depend sensitively on the initial 
conditions. In the low amplitude regime (below a particular critical amplitude 
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(a)' (b)' 

Fig. 3. Particle trajectories (main figures) and the corresponding phase space plots 
(insets) showing different locking behaviour of the particle. For a = 0.5931 (fig. a), 
the particle gets locked in the positive sense while for a = 0.592 (fig. b), the particle 
gets locked in the negative sense, r = 10, x(0) = n/2, v(0) = for both cases. 

depending on the frequency) the particle, after initial transients, gets invari- 
ably trapped to a periodic attractor in some potential minima x = (2n + l)7r/2 
(n = 0,±1,±2,...) and executes periodic motion within the potential well 
(Fig.l). On the average these cases give zero net current as is evident from 
the phase space trajectory (Inset, Fig.l). Only beyond certain amplitude a (for 
a given initial condition and drive frequency u), the particle starts evolving 
giving rise to either periodic or chaotic dynamics. 

Fig. 2 shows the mean velocity < v > for a typical single particle trajectory 
(with the same initial condition (x(0) = tt/2 and v(0) = ) as a function of 
drive amplitude a and period r = 10, where r = 2n/u. There are numerous 
changes in the direction of current as the amplitude a is set at different values. 
For example, whereas a = 0.5931 gives a positive < v >, a — 0.592 for the same 
initial conditions gives mean velocity -< v >. This behaviour can be traced 
to amplitude dependent phase locked periodic behaviour of particle motion 
as shown in Figs.3a,b. As the border between the basins of these attractors 
are crossed with the variation of the amplitude of the drive the single particle 
shows a change in the direction of < v > [2~l~f2l"] . 

Also, the sensitivity of the trajectories to initial conditions can be inferred by 
the following observation. For a = 0.7 and r = 10, the intial condition v(0) = 
0, x(0) = 1.38309 gives a transporting phase locked trajectory in the positive 
direction as in Fig. 3a, whereas v(0) = 0, x(0) = 1.383093 gives locked periodic 
nontransporting trajectory as in Fig.l and v(0) = 0, x(0) = 1.38309704 results 
in a transporting phase locked trajectory in the negative direction as in Fig. 3b. 
The sensitivity of the initial conditions demands that meaningful net particle 
velocity for a given a and u be obtained by averaging over a large number of 
particle trajectories corresponding to various initial conditions. 
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As the nature of the particle trajectory depends critically on the chosen initial 
condition, for realistic averages we calculate the velocity (Eq.6), averaged over 
an ensemble of iV (typically equal to 500) particles. The particles are taken 
with identical initial velocities (i^ n $ = 0), but with initial positions uniformly 
distributed between two consecutive potential maxima [20]. The particle tra- 
jectories, and also the phase space diagrams, are studied to reveal the charac- 
teristics of the particle dynamics. The bifurcation diagrams are obtained by 
recording the values of x(t p ) at times t p = n p r, for each value of the control 
parameter (the amplitude of the external drive), where r = 2ii/u is the period 
of the external drive and n p is an integer. 

3.2 Ratchet current and bifurcation diagrams 

0.8 r 

0.6 - 

0.4 ■ 

0.2 - 

A 
V 

-0.2 - 

-0.4 - 

-0.6 ■ 

0.2 0.3 0.4 0.5 0.6 0.7 0.2 0.3 0.4 0.5 0.6 0.7 

(a) 2 (bf 

Fig. 4. Variation of ensemble averaged velocity ((v)) with amplitude of drive a 
(a). Fig. 4b. shows the corresponding bifurcation diagram. Inset of (b) shows the 
amplitude at which the current reversal occurs in Fig. 4a. r = 10, 7 = 0.12, A = 0.9, 
(f> = 0.35. 

On ensemble averaging, ratchet currents (Fig. 4a) are obtained in the determin- 
istic inertial inhomogeneous system with a symmetric and periodic potential. 
The averaging is done with the same set of parameters as in Fig. 2. The system 
shows current reversals as a function of the amplitude of the drive. However, 
most of the single particle current reversals appearing in Fig. 2 disappear. This 
is due to the fact that different initial conditions select different attractors for 
a particular amplitude, all the other parameters remaining the same. So, at a 
particular amplitude of drive, some particles of the ensemble may move in the 
positive direction and some may move in the negative direction with differ- 
ent multiples of the fundamental locking velocity [21] ■ Others may show a 
chaotic dynamics giving rise to either zero velocity or some non-zero velocity 
on an average. The corresponding bifurcation diagram (Fig.4b) reflects this 
argument. 

Fig.4b shows that the particle exhibits regimes of periodic and chaotic dynam- 




8 



0.03 




(a) 3 (b) S 

Fig. 5. Variation of ensemble averaged velocity ((v)) with amplitude of drive a for 
r = 100 (a) and r = 1000 (b); 7 = 0.12, A = 0.9, (f> = 0.35. 

ICS cLS £L function of the control parameter a. The current reversal at a = 0.6615 
in Fig. 4a is due to a bifurcation from a chaotic to a periodic regime (inset of 
Fig. 4b ) which agrees with Mateos conjencture [20]. However, there could be 
many such bifurcations but without an accompanying current reversal [19]. 

3.3 Frequency dependence of particle dynamics 

The nature of the particle dynamics i.e. the nature of the current and its re- 
versals changes with the driving frequency u = 2-k/t. A comparison between 
the current profile for three periods r = 10 (Fig.4), r = 100 and r = 1000 
(Fig. 5) reflects this fact. For different frequencies, particle current starts de- 
veloping at different amplitudes. This is because the amplitude of the drive 
at which the particle starts evolving out of a single-well trapped-state differs 
with frequency. For very low frequencies of the drive r = 1000 (Fig.5b), par- 
ticle currents are obtained only beyond the critical tilt (a = 1.0 for our case) 
whereas for r = 10 current starts developing essentially around a = 0.55. Also, 
whereas for the slower drive r = 1000 there is no current reversal in the range 
of drive amplitudes shown, for r = 10 and 100 there are several reversals of 
the direction of current in the same range of a. Thus, the mean velocity is 
expected to show a very complex behaviour in the (a — u) space [34J. 

3.4 Role of frictional inhomogeneity and dependence on initial driving phase 

The particle dynamics in the present dissipative system with different degrees 
of frictional inhomogeneity (A), for <fi = 0.35, starting with A = is investi- 
gated. We note that A in the present work plays a role analogous to the po- 
tential asymmetry in the uniform friction case. Surprisingly, at certain drive 
amplitudes, even with zero inhomogeneity, non-zero (steady state, t — > 00) 



9 



0.8 i 1 1 1 1 1 1 3 




%.5 0.55 0.6 0.65 0.7 0.75 "'i).5 0.55 0.6 0.65 0.7 0.75 

(a)" (bf 

Fig. 6. (f) versus a plot for a single particle with A = (a) and the corresponding 
bifurcation diagram (b). r = 10, x(0) = 7r/2, v(0) = 0. 
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Fig. 7. Plot of ((f)) versus a with zero inhomogeniety (A = 0) for initial phase of 
the external drive 4>o = (lower dotted line) and <tj = tt (upper solid line). Middle 
line shows ((f)) after averaging over initial phases. Zeroline is drawn for reference. 
T = 10, A = 0.9, 7 = 0.12, 4> = 0.35. 

current is obtained, both at the single realisation (Fig.6) case, namely (v), 
as well as when averaged over a large number of realizations, ((v)) (Fig. 7). 
Instead, one would have expected the dissipative system to lose the initial 
memory by the time the transients get settled down. The velocities (v) wildly 
change with drive amplitude (Fig. 6(a)) but are finite, ((f)) are not only finite 
but they vary smoothly with drive amplitude. This is obviously an erroneous 
result. Moreover, the bifurcation diagrams with inhomogeneity (Fig. 4(b)) and 
without inhomogeneity (Fig.6 (b)) on comparison appear very similar and do 
not provide any clue to the error. The error could only come from lack of 
averaging over the initial phase 0o of the zero-mean external periodic drive 
F(t). The relevance of averaging over initial phase has earlier been discussed 
for Hamiltonian systems 

In Fig. 7, the ensemble averaged current ((f)) are shown for two values of the 
initial phase O — and tt of the external drive F(t) of period r = 10. Clearly, 
the initial bias due to a fixed value of O determines the steady state current. 
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Fig. 8. Variation of {(v)) with the asymmetry parameter A for initial phase of the 
external drive 4>q = 0(lower dotted line) and cfto = 7r(upper solid line); r = 10, 
a = 0.67, 7 = 0.12, <f> = 0.35. 

In order to get a sensible ratchet current one needs, therefore, to average over 
the initial phase 0o too. The figure also indicates that averaging should be 
carried out over at least two values of 0o which differ by 7r . 

Interestingly, this averaging is required only for small inhomogeneities (Fig. 8). 
For large asymmetries {(v)) does not seem to depend on the initial phase O - 
It turns out that the critical inhomogeneity A = A c ~ 0.6 (Fig. 8). In the case 
of uniform friction the critical potential-asymmetry b ~ 0.25 which is what 
has been used incidentally by Mateos[20j and others corresponding to their 
parameter values. As the value of A is changed, the nature of the particle 
trajectories changes qualitatively. 

In Fig.9 we present the typical trajectories x(t) along with the stroboscopic 
phase points at regular time intervals t = nr,n = 1,2,3, .. for various values 
of the inhomogeneity parameter A for amplitude a = 0.67, period r = 10 and 
initial phase 0o — of the external field F(t). For A = 0.9 all the trajectories 
are chaotic but nontransporting giving zero ratchet current. For A = 0.8 the 
trajectories are chaotic consisting of longer quasiperiodic transporting trajec- 
tories in the positive direction but interrupted by shorter similar trajectories 
in the negative direction giving finite ratchet current in the positive direction. 
For A = 0.7 the current {(v)) is close to the maximum (=2tt/t) corresponding 
to the fundamental velocity. In this case all the trajectories are quasiperiodic 
transporting in the same positive direction. The corresponding stroboscopic 
plot shows their chaotic nature. As the A is lowered from 0.7 to 0.6 we obtain 
a mixture of transporting quasiperiodic trajectories and transporting regular 
periodic trajectories. However, as the value A = 0.6 is reached their directions 
of average motion gradually become equal giving zero current at A = 0.6. As 
mentioned earlier, ((f)) is independent of 0o for the considered parameters. 
For values of A < 0.5 all the trajectories are transporting regular periodic 
with same average slope but differing in direction. For 0o — there are more 
negative-slope trajectories and vice-versa for O — n - The five-paneled Fig.9 
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Fig. 9. Particle trajectories (a), and the corresponding stroboscopic plots (b), for 
different values of inhomogeneity parameter A (A = 0.9 (top) to A = 0.5 (bottom), 
r = 10, a = 0.67, 7 = 0.12, (j) = 0.35. Trajectories are shown for ten different sample 
initial conditions. Though the trajectories appear to be straight line in the scale of 
the figures, they are oscillatory as in Fig. 3 

gives a qualitative mechanism for the explanation of Fig. 8 and underlines the 
importance of averaging over the initial phase <po at smaller inhomogeneities. 

3. 5 Role of noise and low amplitudes of drives 

In deterministic ratchets, because of the complex structure of the basins of at- 
traction and the presence of multiple attractors in phase space [38], the nature 
of the current is sensitively dependent on the initial conditions. The presence 
of noise in the system, however, makes the particle dynamics independent of 
the initial conditions allowing the particle to explore all the existing attrac- 
tors. As a result the chaotic behaviour of the system disappears. The average 
velocity shows a relatively smoother variation with the amplitude of exter- 
nal drive than in the deterministic case. In Fig. 10, the variation of v with a 
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Fig. 10. Plot of v versus a in the presence of noise, for three values of r; T = 0.4, 
7 = 0.12, A = 0.9, <j) = 0.35. 

is shown for three different values of drive frequencies, keeping all the other 
parameters fixed. 

With noise, appreciable ratchet current is obtained even at low amplitudes, 
where there was no current in the deterministic case (Fig. 10). This is be- 
cause, the presence of noise aids the particles to overcome the potential barrier 
which it could not have done otherwise. Hence the presence of noise plays a 
non-negligible role in obtaining ratchet current in the low amplitude regime, 
contrary to what has been mentioned by others [20]. The same effect of noise 
has been shown in a different system using the matrix continued fraction 
method [41 J. 



4 Conclusion 

In this work the deterministic ratchet current was obtained in a periodically 
driven symmetric potential in the presence of frictional inhomogeniety (most 
of the earlier works were with a driven asymmetric potential). The particle 
dynamics shows similar characteristic as in an asymmetric potential. Initial 
conditions and the control parameters play a major role in shaping the nature 
of the dynamics. For low amplitudes of drive the initial phase of the drive 
is found to play a non-negligible role. At these amplitudes, we show that 
averaging over initial conditions and initial phases of the external drive is 
essential to obtain realistic averages. With noise, ratchet current is obtained 
at amplitude regimes where there is no current in the deterministic systems. 
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